Magnetic anisotropy in Li-phosphates and origin of magnetoelectricity in LiNiP04 
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Li-based phosphates are paradigmatic materials for magnetoelectricity. By means of first- 
principles calculations, we elucidate the microscopic origin of spin anisotropy and of magnetoelectric 
effects in LiNiP04. The comparison with LiCoP04 reveals that Co-d^ and Ni-d^ electronic clouds 
show distinct orbital shapes, which in turn result in an opposite trend of the local spin anisotropy 
with respect to the surrounding Oe cages. Due to magnetic anisotropy, the Ni-based phosphate 
shows a peculiar "angled-cross" spin ground-state, which is responsible for magnetoelectricity. In 
this respect, we show that, under a magnetic field Hx, an electronic polarization Pz arises, with an 
estimated linear magneto-electric coefficient in good agreement with experiments. 

PACS numbers: Vahd PACS appear here 



Olivine phosphates LiCoP04 and LiNiP04 are at- 
tracting large interests, due to their peculiar magne- 
toelectric (ME) effect [i.e. the control of ferroelec- 
tric (magnetic) properties via a magnetic (electric) field] 
as well as their application for electrodes in recharge- 
able Li batteries [J [2[. Recently, ferrotoroidic do- 
mains and antiferromagnetic domains have been indepen- 
dently observed in LiCoP04 by using second harmonic 
generation. [3] The toroidal moment, generated by a vor- 
tex of magnetic moments, is considered as source of a 
novel ferroic order, closely related to ME effects [1, [H as 
well as to multiferroicity^Gj (^-e. coexistence of long-range 
magnetic and dipolar orders). In fact, LiCoP04 shows a 
nonzero linear ME coefficient, axy and ayx, at low tem- 
perature, consistent with the toroidal moment T nearly 
parallel to z axisQ, whereas LiNiP04 shows axz and 
<^;zcc0- Although magnetoelectricity has been macro- 
scopically investigated by means of Landau theoryjl, 
its microscopic origin has not been clarified yet nor first 
principles calculations aimed at investigating ME effects 
in phosphates exist in the literature. In LiCoP04, anti- 
ferromagnetic (AFM) spins were found to be along the b 
axis[10] or uniformly rotated from this axis by 4.6 
In LiNiP04, it was proposed the collinear AFM spins 
to lie along the c axis but recently a non-collinear 
structure in a "spin-cross" style [1, [l2[ was suggested, in 
order to explain the butterfly shape of the ME hysthere- 
sis curve |13|. 

In this letter, within Density-functional theory (DFT) 
we investigate the magnetic anisotropy for Li-phosphates 
LiTMP04 (TM = Ni, Co) and we are able to correlate the 
orbital degree of freedom with the calculated anisotropy 
and spin-configuration. The discussion on magnetic 
anisotropy is particularly important for LiNiP04, where 
we find, as magnetic ground-state, a peculiar spin-cross 
configuration that leads to magnetoelectricity; finally, we 
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clarify the microscopic origin of ME effects, our calcu- 
lated ME coefficient being in quantitative agreement with 
experiments. 

Methodology and structural details. DFT simulations 
were performed u sing the VASP code [3] and the PAW 
pseudopotentials [l5[ within the GGA+t/ formalism p!6j 
{U=5 eV and J=0 eV for Ni d-states). Other values of 
U=2 and 8 eV were also tested. The cut-off energy for 
the plane-wave expansion of the wave-functions was set 
to 400 eV and a k-point shell of (2, 4, 4) was used for the 
Brillouin zone integration. Lattice parameters were fixed 
as experimentally observed [13, US] • The internal atomic 
coordinates were fully optimized in a (fully-compensated) 
AFM configuration, keeping Si = S2 = —Ss = —S4 
configuration. In the orthorhombic Prima (DIq) para- 
magnetic space group, these four Co/Ni sites are re- 
lated by eight symmetry operations; all the rotations 
and mirror reflections accompany translations, so that 
the Co/Ni sites are slightly deviated from high-symmetry 
positions. Optimized coordinates of the four Co/Ni 
ions are: Col/Nil(l/4+e, 1/4, -S), Co2/Ni2(3/4+e, 1/4, 
1/2+J), Co3/Ni3(3/4-e, 3/4, 6), Co4/Ni4(l/4-e, 3/4, 
1/2-5), where e=0.0268 and 5=0.0207 for Co, e=0.0249 
and 5=0.0152 for Ni. The spin-orbit coupling (SOC) 
term was computed self-consistently inside each atomic 
sphere, with a radius of 1 A. The electronic polarization 
P was calculated using the Berry phase method [Tol]. 

Magnetic anisotropy: LiCoPO^ vs LiNiPO^. As 
shown in FigiH each Co ion is surrounded by highly dis- 
torted oxygen octahedra, so that the partially filled t2g 
shell shows a dumbbell-shaped charge distribution, axi- 
ally elongated along the iii direction (along ml-m2 0-Co- 
O bonds). This peculiar shape is expected to induce a 
strong local magnetic anisotropy via the SOC term. On 
the other hand, in LiNiP04, the Ni-(i^ orbital shows a 
more isotropic sphere-like shape. 

Both the global and local magnetic anisotropy were 
investigated by rotating four spins simultaneously, keep- 
ing a collinear AFM coupling. The global magnetic 
anisotropic energy (MAE) was calculated by differ- 
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FIG. 1: (a) Isosurface of charge density of minority-spin Co- 
t2g states (within an energy range up to 1 eV below Ef)- Co- 
O bond lenghts (A) are reported and denoted as /,ml,m2,s. 
Easy planes for Co local spins and the corresponding normal 
vectors (hard axis) iii are also shown at Col and Co2 sites, 
(b) Section of charge density of Co t2g electrons and magnetic 
easy and hard axis in the ac plane in LiCoP04. (c) Same as 
(b) for LiNiP04. Unfilled large arrows show hard (rii) and 
easy (n^) axis in ac plane for Co and Ni spins respectively. 



ences in the total energy with different spin-orientations, 
whereas the local anisotropic energy was evaluated as 
proportional to the expectation value of the SOC en- 
ergy: £;soc = (^^^^s)[20] (integrated in each atomic 
sphere) . 

Fig. [2] (a) and (b) shows the MAE of Co spins in 
the ac and ab planes, respectively. The global easy axis 
is the b direction, in agreement with the experimental 
suggestion that spins should be aligned along the b axis 
(possibly with a slight rotation). The MAE is rather high 
(more than 10 times larger than the orbitally-ordered 
Mn spins in TbMn03[21|]). As shown by the spin angle 
dependence of Esoc (cfr Fig. [2] (a)), the easy/hard axial 
direction depends on each Co site. Etot can be fitted by 
a conventional quadratic anisotropic term for Si spins: 
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FIG. 2: Total energy AEtotai (black solid line) and AEsoc 
(dashed colored lines) at each Co site vs the direction of 
collinear Co spins (a) in ac plane and (b) in ab plane. Cor- 
responding curves for Ni spins are shown in (c) and (d) for 
the ac and ab planes, respectively. Note the different energy 
scales of panels (a) and (b) compared to (c) and (d). The 
energy with spins along the a axis is taken as reference. 



I)(Si/|Si| -iii)^, where iii is the site-dependent hard axis, 
described as ni = ns = (cos a, 0, sin a), n2 = n4 = 
(cos a, 0, — sin a) (see. Fig(T]). According to the fitting, 
D=7.39 meV and a=35.48°: this implies the hard axis 
iii to be nearly along the ml bond, so that Co spins are 
aligned in an easy plane perpendicular to iii. In terms of 
site-dependent anisotropy, we note that the stable AFM 
spin ordering Si = S2 = —Ss = — 5*4 is different from 
the orbital ordering Li = L^^ L2 = L4 which causes a 
local anisotropy. In such a configuration, all four spins 
are allowed to lie in the easy plane only when collinear 
AFM spins are pointing along the b axis. Therefore the b 
axis is the easy axis in the collinear AFM configurations 
for Co spins. 

Our calculations show for Ni spins an opposite trend 
of EsoCi with respect to the Co spin, both in the ac 
and ab planes (cfr Fig. [2|); moreover, the MAE is al- 
most two orders of magnitude smaller, consistent with 
the spherical electronic cloud. Here the b axis is the 
global hard axis. Assuming the easy axis to lie in the ac 
plane and by fitting Etot to the anisotropic term, we es- 
timate D=-0.50 meV and a=46.7° Here a is close to 7r/4 
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FIG. 3: Total energy AEtotai (black solid line) and AEsoc 
(dashed red line) at each Ni site vs the direction of non- 
collinear Ni spins in the "spin-cross" configuration (shown in 
the inset). The minimum is marked by a vertical blue arrow. 
The dashed line is a function proportional to 1 — cos^, as a 
guide to the eye, to outline the "asymmetrical" behaviour of 
the DFT data. 



so that two opposite contributions from local anisotropy 
nearly cancel out (Fig. [21(c)) and induce a small in-plane 
anisotropy. The negative value of D implies the existence 
of an easy axis iij of Ni spins in the ac plane (see. FigHj), 
n'^ = rig = (— sina, 0, cosa), 112 = 114 = (sin a, 0, cos a) 
which is close to the iii direction at Co spins. Therefore, 
the c axis is the global easy axis for Ni spins if a collinear 
AFM configuration is assumed. However, in order to sta- 
bilize the local anisotropic term, the spins are expected 
to slightly tilt with respect to the c axis. This deviation, 
denoted as "angled cross" spin configuration, has been al- 
ready discussed by Chupis[9| in terms of Landau theory 
and suggested to play an important role for magnetoelec- 
tricity. This issue has been carefully investigated here by 
tilting the spins by an angle from the c axis. Indeed, as 
shown in Fig. [3l the non-collinear spin structure in the ac 
plane with ^ ~ 1° gives the lowest energy. We considered 
a generic Hamiltonian for Si (i = 1 to 4) spins, including 
a Heisenberg term, an relativistic anisotropy term and 
a Zeeman-like term (see next paragraph, where a finite 
external H-field will be introduced): 
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In this H=0 case, a delicate balance of the first two terms 
occurs: As the spin is tilted towards the local easy axis, 
the SOC-term is stabilized, whereas the Jij coupling be- 
comes unstable. As a result, the equilibrium occurs in a 
non-collinear "angled-cross" spin-configuration. 

Magnet oelectricity in LiNiPO^ . The existence of a lo- 
cal magnetic anisotropy is particularly important in the 
context of magnetoelectricity. In order to investigate ME 
effects in LiNiP04, hereafter we focus on the change in 
the spin configuration under a magnetic field with re- 
spect to the "angled cross" ground-state. The exchange 
coupling constants Jij were evaluated by total energy 
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FIG. 4: (a) Spin configuration in the ac plane under an ap- 
plied field Hx . (b) Magnetic- field induced electric polarization 
along z (with and without SOC, as well as considering D=0 
in Eq.(l), see text) in LiNiP04. The magnetization Mx is 
also shown (left-trianges, referred to the right 2/- axis) Inset: 
spin angle 01 and 02 (in degrees) vs Hx (T) 



differences fitted to Eq. ([T]), considering several AFM 
configurations. Our calculated values are J12 = J34 = 
-O.llSmeF, Ji3 J24 lAGmeV, J14 = J23 3.92meV 
(the positive sign means AFM coupling). The dominant 
Ji4 keeps the spins in the AFM configuration, whereas 
J12 and Ji3 are responsible for removing the degeneracy 
of some AFM configurations [8j, i.e. the AFM spin config- 
uration 81=82=— Ss=—S4 gives the most stable energy. 
The J values are consistent with the fact that four Ni 
sites are separated into two pairs (Nil, Ni4) and (Ni2, 
Ni3) by a Li intercalated layer, so that different pairs are 
expected to be weakly coupled. The anisotropy coeffi- 
cient D is also obtained from DFT, as reported above. A 
nonzero applied magnetic field, Hx^ is expected to tilt the 
Ni spins from the angle-crossed spin configuration, with 
angles from the c axis, denoted as Oi and ^2 -[HI] (Here 
we assume the relation Oi = O4 and O2 = Os by consid- 
ering the symmetry under the magnetic field.) By using 
DFT parameters, Eq.(l) is easily minimized by Newton 
method with Oi and O2 angles varying by few degrees in 
a way proportional to Hx (cfr FigHlb), inset). 

Constraining the spin angles as obtained from the New- 
ton minimization for each H-field, we evaluate the elec- 
tronic polarization. A small but finite Pz is obtained, 
confirming magnetoelectricity. As shown in Fig. [H all 
spins have positive x components for Hx >5T; however, 
the sign of the moment doesn't affect the induced P^, 
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which is rather hnearly proportional to H^^ consistent 
with the experimentahy observed behaviour at low tem- 
perature. From the slope of the P-H curve from H=0 
to 10 T, the linear ME coefficient azx is estimated as 
0.59 iiC/im?T — {).lAps/m^ comparable with the low- 
temperature experimental value = 1.5p5/m.[l2| 

Let us now investigate the microscopic origin of mag- 
netoelectricity. Electric dipoles are found to be mainly 
induced by symmetric exchange [23|. often referred to as 
"inverse" Goodenough-Kanamori (iGK) interaction [23] 
in Nil-0-Ni4 (or Ni2-0-Ni3) bonds. For smah H, the 
other conventional mechanism for magnetically-induced 
polarization - related to antisymmetric exchange and of- 
ten labeled as "inverse Dzyaloshinskii-Moriya (iDM)" - 
is negligible. However, for high magnetic fields, the iDM 
becomes sizeable (cfr FigHlb), where the comparison be- 
tween Pz values including - or not - the SOC term is 
reported). In this context, we note that the existence of 
a local anisotropic term is important to induce P\ As- 
suming I) = in Eq.(l), so that 01 = 02^ we find that the 
induced Pz is less than half with respect to the nonzero 
D case (cfr Figglb). 

Finally, we estimate toroidal moments based on Ref. ^ . 
The collinear AFM spin configuration in LiCoP04 in- 
duces (0, 0, -1.52)/iBA , whereas the noncollinear spin 
configuration (at H=OT) in LiNiP04 induces (0.19, 0.94, 



0) /i^A per unit cell. We observe that the spin crossed 
state causes a small Ta^-component; however the Ty com- 
ponent, coming from the compensated AFM spin config- 
uration, is relevant to the ME effect with finite Pz and 

In conclusion, we have presented a careful investigation 
of the spin-anisotropy in Li-based transition-metal phos- 
phates. We found that LiNiP04 shows a peculiar ground- 
state spin configuration: due to a site-dependent mag- 
netic anisotropy, it shows non-collinear spins arranged 
as "angled-cross"-like. This is relevant in the context of 
magnetoelectricity: due to the ground-state noncollinear- 
ity, the spin configuration induced by an applied mag- 
netic field leads to a net electric polarization, as shown 
by a realistic first-principles estimate of the magnetoelec- 
tric coefficient in agreement with experiments. Our re- 
sults suggest a possible avenue for elucidating the origin 
of the ME effects in other compounds. 

Acknowledgments. We thank Claude Ederer for helpful 
discussions. The research leading to these results has re- 
ceived funding from the European Research Council un- 
der the EU Seventh Framework Programme (FP 7/2007- 
2013) / ERC grant agreement n. 203523. Computational 
support from Caspur Supercomputing Center (Rome) 
and Cineca Supercomputing Center (Bologna) is grate- 
fully acknowledged. 



[1] Amine, K; Yasuda, H; Yamachi, M, Electrochemical and 

Solid State Letters 3, 178 (2000). 
[2] Le Bacq, O; Pasturel, A; Bengone, O, Phys. Rev. B, 69, 

24 (2004). 

[3] B. B. Van Aken, J. P. Rivera, H. Schmid and M. Fiebig, 

Nature 449 702 (2008). 
[4] Schmid, H. Ferroelectrics 252, 41750 (2001). 
[5] C. Ederer and N. A. Spaldin, Phys. Rev. B 76 214404 

(2007). 

[6] N. A. Spaldin and M. Fiebig, Science 309, 391 (2005). 
[7] Rev. Gen. Electr. 80 143 (1971). 
[8] D. Vakinin et. al., Phys. Rev. B 65 224414 (2002). 
[9] I. E. Chupis, Low Temperature Physics 26 419 (2000). 
[10] R. P. Santoro, R. E. Newnham and S. Nomura, J. Phys. 

Chem. Solids 27 655 (1966). 
[11] R. P. Santoro, D. J. Segal, and R. E. Newman, J. Phys. 

Chem. Solids 27 1192 (1966). 
[12] T. B. S. Jensen, et. al., Phys. Rev. B 79, 092412 (2009). 
[13] J.-P Rivera, Ferroelectrics 161, 147 (1994). 
[14] G.Kresse and J.Furthmiiller, Phys.Rev.B 54, 11169 
(1996). 

[15] P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
[16] V.LAnisimov, F.Aryasetiawan and A.LLichtenstein, J. 
Phys.: Cond. Mat. 9, 767 (1997). 



[17] D. Vaknin,l J. L. Zarestky, L. L. Miller, J.-P. Rivera and 
H. Schmid, Phys. Rev. B 65, 224414 (2002). 

[18] L Abrahams and K. S. Easson, Acta Cryst. C49, 925- 
926. (1993). 

[19] R.D.King-Smith and D.Vanderbilt, Phys. Rev. B 47, 
1651 (1993); R. Resta, Rev. Mod. Phys 66, 899 (1994). 

[20] J. Kiibler, Theory of Itinerant Electron Magnetism (Ox- 
ford Science Publications, New York, 2000). 

[21] H. J. Xiang, Su-Huai Wei, M.-H. Whangbo and Juarez 
L. F. Da Silva,Phys. Rev. Lett. 101, 037209 (2008). 

[22] The effect of an applied magnetic field is taken into ac- 
count only through local spin directions {i.e. neglecting 
dynamical effects by the field). 

[23] S. Picozzi, L A. Sergienko, K. Yamauchi, B. Sanyal and 
E. Dagotto, Phys. Rev. Lett. 77, 227201 (2007); S. Pi- 
cozzi and C. Ederer, J. Phys. Cond. Matt. 21, 303201 
(2009). 

[24] K. Yamauchi, S. Picozzi, J. Phys. : Condens. Matter 21, 
064203 (2009). 

[25] H.Katsura, N.Nagaosa, A.V.Balatsky, Phys. Rev. Lett. 
95, 057205 (2005); LA. Sergienko et al. and E.Dagotto, 
Phys. Rev. B 73, 094434 (2006); M.Mostovoy, Phys. Rev. 
Lett. 96, 067601 (2006). 



